The QCD chiral transition temperature in a Dyson-Schwinger-equation context 



M. BlanlQ and A. Krassnigjt] 

Institut fur Physik, Karl-Franzens-Universitat Graz, A-8010 Graz, Austria 
(Dated: April 30, 2010) 

We analyze the chiral phase transition with the help of the QCD gap equation. Various models 
for the effective interaction in rainbow truncation are contrasted with regard to the resulting chiral 
transition temperatures. In particular, we investigate possible systematic relations of the details 
of the effective interaction and the value of T c . In addition, we quantify changes to the transition 
temperature beyond the rainbow truncation. 
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I. INTRODUCTION 

The structure of the phase diagram of quantum chro- 
modynamics (QCD) has recently received a lot of atten- 
tion. Theoretical progress in various directions is well 
in line with past, present, and upcoming experimental 
programs e.g. at CERN SPS, RHIC, LHC, and FAIR. 
Among the key features of the phase diagram is the chiral 
phase transition at nonzero temperature and vanishing 
chemical potential. This transition has been and is stud- 
ied extensively by various different approaches with an 
emphasis on the critical behavior of the order parameter 
as well as on the actual value of the transition tempera- 
ture T c . Lattice QCD calculations for example over the 
past years yielded T c at temperatures between 145 and 
200 MeV with a recent tendency towards a commonly 
found value [l]-0] • In various model approaches in mean- 
field approximation and beyond one finds a range from 
90 to 160 MeV, see e.g., [§[ and references therein. 

Although an understanding of the universal properties 
of QCD is of great importance, it is also rewarding to 
investigate the value of T c which is not universal and 
can therefore be used to draw conclusions about specific 
details of the underlying dynamics of quarks and gluons. 
More precisely, specific models and their assumptions can 
be falsified in the approach presented here by a calcula- 
tion of T c in QCD. 



II. FINITE-TEMPERATURE APPLICATIONS 
OF DYSON-SCHWINGER EQUATIONS 

A nonperturbative continuum approach to QCD is pro- 
vided by the Dyson-Schwinger equations (DSEs) (9l.fTfjj. 
which can be applied straight-forwardly also for finite 
temperature and chemical potential [ll|. The DSEs are 
well-suited to study the properties of QCD's elementary 
degrees of freedom, quarks and gluons. In this frame- 
work, various aspects of chiral symmetry and its dynam- 



ical breaking have been studied some time ago via the 
DSE of the quark propagator (also referred to as the 
"QCD gap equation") at zero temperature (see, e.g. [II] 
and references therein) and in connection with the phase 
diagram of QCD (see, e.g. [Tl|] and references therein). 
The same framework also allows the consistent and co- 
variant study of bound-state, i. e. meson and baryon, 
properties. At zero temperature, a level of sophistica- 
tion beyond simple QCD modeling has been reached (see 
e.g. [L3l - ll5j ] and references therein). At nonzero tem- 
perature and density, analogous investigations are com- 
plicated by the additional variables introduced in the 
Matsubara formalism. Nonetheless, investigations of me- 
son masses at and beyond the phase transition at zero 
chemical potential 
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I- 22 1 as well as certain aspects of 
U28I ] have been conducted at vari- 
ous but mostly simpler degrees of sophistication. In the 
present context, our focus is on DSE studies of the chiral 
phase transition in QCD at finite temperature and zero 
chemical potential. 

It was shown some time ago that certain classes of 
models in the QCD gap equation in rainbow trunca- 
tion yield a second-order phase transition with mean- field 
critical exponents [2t| [30] ■ While the main emphasis of 
that investigation was the critical behavior and universal- 
ity class of the phase transition, it also provided the cor- 
responding values for T Cl which ranged between 120 and 
174 MeV among the model interactions investigated. In 
somewhat different setups, using a separable interaction 
kernel in the gap equation [r7l 0, Ho, [Hj] or neglecting 
retardation effects [2lj . one obtaines values in the range 
between 110 and 146 MeV. 

The pointwise behavior of all these interactions as a 
function of the gluon momentum is rather different, al- 
though in all cases it was determined by adjusting the rel- 
evant parameters to meson phenomenology at zero tem- 
perature. Therefore, it is interesting to ask whether there 
is any simple relationship between the momentum depen- 
dence of the interaction and the value of T c . Such a re- 
lationship could e.g. be analogous to that found recently 
in studies of meson properties using a particular form 
of model interaction with a one-parameter setup. There 
the value of the free parameter determines an effective 
range of the intermediate- and low-momentum parts of 
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the interaction. While ground-state properties of pseu- 
doscalar and vector mesons were unaffected by variations 
in the model parameter, masses of excitations of any kind 
(orbital or radial) showed a strong and systematic depen- 
dence on the model parameter [141 132-35] thus identify- 
ing excited-state properties as prime targets to study the 
particular details of the effective interaction in the DSE 
formalism, in particular in the nonperturbative regime. 
In an analogous fashion, the goal of the present work is 
to identify T c as a quantity with the same capability. In 
particular, one can expect a calculation of T c to impose 
even stronger restrictions on model parameters, beyond 
what is at hand at zero temperature. 

The paper is organized as follows: Sec. IIIII briefly 
sketches the quark DSE in QCD at finite temperature, 
Sec. IIVI compares the model interactions used herein, 
Sec. |V] outlines how we compute T c and the critical expo- 
nents. The results are presented in Sees. IVII and lVIll the 
latter focusing on effects from corrections beyond rain- 
bow truncation. Conclusions and an outlook are pre- 
sented in Sec. IVIIII All calculations are performed in 
Euclidean space. 



III. QCD GAP EQUATION AT FINITE 
TEMPERATURE 



We use the quark DSE to study the chiral phase tran- 
sition of QCD. It is an inhomogeneous nonlinear integral 
equation for the dressed-quark propagator, whose inverse 
at T = is of the form 



S- 1 ( P ) = t 1 -pA(p 2 ) + B(p 2 ), 



(1) 



where the four-vector p denotes the quark momentum 
and A and B are scalar functions of the momentum- 
squared. At nonzero temperature the presence of the 
heat bath implies an additional four vector, which leads 
to more structures in (JI]) [3a|. In the Matsubara formal- 
ism one arrives at 

S , ~ 1 (p,w fe ) = ij-p A(p 2 ,Lo k ) +B(p 2 ,uj k ) 

+ z7 4 w fe C(p 2 ,w fe ) + 7 • p -f4^kD(p 2 ,waJ2) 

Now the four scalar functions (A, B, C, D) depend on the 
quark 3-momentum squared p 2 and the quark's fermion 
Matsubara frequency u> k = 2irT(k + 1/2), k 6 Z. 
The QCD gap equation in pictorial form reads 



K^r 1 = (— 

S(p,KJ k ) So(p>'fc) 



) ~ 1 + | , r " (?? ' Uk ' ^ " ;) ' ^ 

S(<f,W() 



where So(p, uik) is the free quark propagator, D^ v {p — 
q, LUk — u!i ) the renormalized dressed gluon propagator 
and r„(p, Wfe; q 7 ui{) the renormalized dressed quark- gluon 
vertex. Solving for S(p,uik) one needs expressions for 
the quark-gluon vertex as well as the gluon propagator. 



TABLE I: Parameter sets D and oj used for the different in- 
teractions MN Ep, Eq. {TDJ; AWW [3^|. Eg. (fn): MT pj. 
Eq. (Tl2)) ; MR [4(]| ]. Eq. (fl3|) : together with corresponding ob- 
servables in T = meson studies, where available. Numbers 
are in GeV except for D (GeV 2 ) and the chiral condensate, 
whose dimension is given explicitly. 



Model 




D 


-{qq)o 


m n 


m e 


MN 




U.JUIO 


('lie: n P V~l 3 


0.14 


n 77 


AWW1 


0.3 


1.47 


(.245 GeV) 3 


0.135 


0.745 


AWW2 


0.4 


1.152 


(.246 GeV) 3 


0.135 


0.748 


AWW3 


0.5 


1.0 


(.251 GeV) 3 


0.137 


0.758 


MT1 


0.3 


1.24 


(.243 GeV) 3 


0.139 


0.747 


MT2 


0.4 


0.93 


(.242 GeV) 3 


0.139 


0.743 


MT3 


0.5 


0.744 


(.239 GeV) 3 


0.141 


0.724 


MR1 


0.3 


0.78 


(.241 GeV) 3 


0.139 




MR2 


0.4 


0.78 


(.250 GeV) 3 


0.139 




MR3 


0.5 


0.78 


(.255 GeV) 3 


0.139 





These can in principle be obtained from their own DSEs 
self-consistently, which is clearly beyond the scope of the 
present study. Instead we use various model Ansatze in 
a rainbow truncation of the gap equation, as specified 
below. Defining the gluon four momentum as (fc, fi) := 
(p— q,uJk — ^i), the second term on the right-hand side 
of (|3]), the quark self energy E(p, ujk), is given by 



T 



' l,q 4g 2 D^(k,n)x 



(2tt) 3 3' 



j v S (q, uji )r M {p, u) k ;q,u>i). 



(4) 
(5) 



The gap equation is subject to the renormalization con- 
dition that for a renormalization scale ( and p 2 +uJq = Q 2 

S^ 1 {p,uj ) =ij -p + ij4,ujQ + m(C). (6) 

Note that this condition does not contain a term 7 • 
pjiLVk, i.e., the scalar function D(p 2 ,u>k) will be power- 
law suppresed in the ultra-violet jllT |. 

In the following, unless otherwise noted, we 
use the rainbow-truncated gap equation by setting 
r M (p, uik ; <f, ) — jp. The finite- T expression for the 
renormalized dressed gluon propagator in Landau gauge 
reads [H ELI 



g 2 D lu/ (k,Q) 



p^(k,n) g(k,n ;mg ) 



0, [i and/or v = 4 

"75*, V,v = i,j = 1,2,3 



Sfj,v - 



k 2 
kyk v 

k 2 + n 2 



p T 



(J) 

(8) 
(9) 



where m g is a Debye mass [42J and Q is an effective in- 
teraction detailed below. 
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IV. INTERACTION MODELS 

At zero temperature covariant bound-state equations 
have been used to study the properties of hadrons in 
various particular setups in order to gain insight regard- 
ing the properties of the underlying fundamental Green 
functions of QCD from a phenomenological point of 
view. More concretely, for the case of mesons the Bethe- 
Salpeter-Equation in ladder truncation together with the 
gap equation in rainbow truncation has been employed 
using many different Ansatze for the quark-gluon inter- 
action to investigate various meson properties, see [lH 
and references therein. In that process, different effec- 
tive interactions have been tested by employing model- 
independent constraints for the interaction and fixing the 
remaining (in our case either one or two) parameters and 
a current-quark mass to phenomenological and experi- 
mental numbers such as the chiral condensate, the pion 
mass, and the pion decay constant. 

In the present work we investigate four different func- 
tional forms of the interaction, referred to in the follow- 
ing by MN (Munczek-Nemirovsky [13), AWW (Alkofer- 
Watson-Weigel 38]), MT (Maris-Tandy [11]), and MR 
(Maris- Roberts [4(|), in order of complexity. The 
zero-temperature effective interaction can be straightfor- 
wardly used at finite temperature by insertion in Eq. (J7J . 
Via s:=k 2 + n 2 



£(s) 



Q(k,fl;m g ) 
MN : 

AWW : 

MT : 

MR : 



k-l.O 



4tt 2 
D — —se 

L0 6 



■s/u 



D 



An' 



(10) 

(11) 
(12) 



Fuv(s) ■■ 



7 OJ" 

+T uv (s), 

4-7T 7 m 7T Fjs) 

l/21n[r+(l + S /A| CD )2] 



(13) 



As given in 
0.5 GeV, r = 
7m = 12/(33 



1, = [1 - exp(-s/[4mf])]/s, m f = 

2 - 1, N f = 4, Aq^ 4 = 0.234 GeV, and 
2Nf). At T = all of these Ansatze pro- 
vide the correct amount of dynamical chiral symmetry 
breaking as well as quark confinement via the absence of 
a Lehmann representation for the dressed quark propa- 
gator. In addition, MT and MR produce the correct per- 
turbative limit of the QCD running coupling, i.e. they 
preserve the one-loop renormalization-group behavior of 
QCD for solutions of the quark DSE. Thus, for MT and 
MR, a renormalization procedure is needed in the gap 
equation, which we implement following Refs. (l8l . l40j |. 
where the necessary details can be found. 

Here, for the interactions in Eqs. (|TT|) . (|T2"j) . and (TT3"1) . 
we investigate different parameter sets (i.e., values of D 
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FIG. 1: (Color online) Plot of the different forms for the ef- 
fective interaction used in this work as simple functions of a 
momentum squared s. The corresponding parameters can be 
found in Tab. [I] Note the different behavior of the interactions 
on all momentum ranges. In particular, the reader should 
keep in mind that the 5-function present in MR, Eq. (|13|) is 
invisible in the double-logarithmic presentation as well as the 
MN interaction, Eq. JT0]|, which is zero for a ^ 0. 



and uj) in analogy to the corresponding bound-state stud- 
ies at T = 0. Table fl] lists these sets of parameters 
together with the respective results for relevant T = 
observables. Once D and w have been fitted to the chi- 
ral condensate, one remains with a u/d current-quark 
mass to be chosen in order to calculate properties of 
light mesons. For each set the light quark mass has 
been adjusted to fit m w ; m g is listed for a first corre- 
sponding result without fixing any further parameters. 
Since our study will remain in the chiral limit, the meson 
masses are included only to complete the justification of 
the model parameters. 

For easy reference and an instructive comparison, the 
effective interaction for all parameter sets is plotted in 
Fig. [TJ Differences are clearly visible in the UV regime. 
In the IR one has to keep in mind that the ^-function 
present in the MR interaction, Eq. (fT3")) is not shown 
here due to the double-logarithmic presentation as well 
as the MN interaction, Eq. (fTU|). which is zero for s / 0. 
The intermediate momentum region is dominated by the 
Gaussian peak whose position and width are character- 
ized by the parameter u> and whose overall strength is 
D. 

The reader may ask whether there is any relation be- 
tween D and u> that could simplify the parameterization. 
Indeed, inspection of Tab. U shows that for the MT in- 
teraction the product u>D is constant over the range of 
values for uj used here, when one aims at a given value 
of the chiral condensate. In addition, without changing 
the value of the u/d current-quark mass, also (and 
approximately also m g ) are unchanged on this domain, 
which can be used to motivate a redefinition of the pa- 
rameters [43j]. As already mentioned in Sec. HH one can 
interpret this pattern with regard to u> such that this 
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FIG. 2: (Color online) Plot of the order parameter Bo vs. 
temperature T, shown exemplarily for model MT2. 



parameter represents an inverse effective range of the in- 
teraction, a dependence on which should be less promi- 
nent in ground-state properties than for excited states. 
However, for the other interaction types, such a simple 
relation among the parameters does not exist and so we 
keep the interactions in their present form. 



V. CHIRAL PHASE TRANSITION 

With all ingredients of the QCD gap equation at hand 
one can obtain a solution numerically and study the 
chiral phase transition temperature T c as well as the 
nature of the transition. The corresponding order pa- 
rameter is the chiral condensate. Equivalently, we use 
B := B(0, wo), since it is easier to access and can be cal- 
culated more accurately. A typical behavior of the order 
parameter for a chiral-limit solution is shown in Fig. [5] 
and indicates a second order phase transition. Indeed as 
mentioned earlier, it has been found previously that the 
rainbow-truncated quark DSE generally yields a second 
order chiral phase transition with mean-field critical ex- 
ponents [ID, [13] ■ We have confirmed this behavior for 
all interactions considered here and thus enlarged the set 
of interaction types it had been shown for by MT and 
AWW. The following account of our procedure includes 
the discussion of critical exponents mainly for complete- 
ness, but also, since as described below, T c can be ex- 
tracted more accurately this way. In the case of a sec- 
ond order phase transition, the order parameter Bo(t, h) 
obeys the scaling laws 



/i=0 



B (t,h) cx (-tf\. 



B (t,h) oc h 1/s 



4 0" 
0+ 



and 



(14) 
(15) 



Here, t 



Sr — 1 is the reduced temperature and h = ^ 
the reduced mass, a measure for the explicit breaking of 
chiral symmetry by a non- vanishing current quark mass 
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FIG. 3: (Color online) Chiral susceptibilities Xt(T) (upper 
panel) and Xh{T) (lower panel) for the MT1 model plotted 
versus temperature T for various values of h. The curves 
between the data points were obtained by interpolation using 
cubic splines. 



to. /? and S are the critical exponents of the phase tran- 
sition. 

Although it is straightforward to use Eqs. (fT4"l) - (fT5")) to 
obtain j3 and 5, this procedure requires to solve the gap 
equation at h — or T — T c . This is numerically difficult, 
such that a direct evaluation does not allow fits to extract 
T c with the necessary precision and reliably observe scal- 
ing. Therefore we exploit further scaling relations and 
use chiral susceptibilities for our analysis. They are de- 
fined by 



XT 



Xh ■= 



-T, 



dB Q (T, h) 



dT 
dB Q (T,h) 



and 



h fixed 



dh 



(16) 
(17) 



T fixed 



The maxima of these quantities for nonvanishing h are 
referred to as pseudocritical points, x V t anc ^ X/T: res P ec_ 
tively. The corresponding pseudocritical temperatures 
are denoted by T|f and Tf^ . They are obtained as the 
maxima of the chiral susceptibilities with respect to tem- 
perature, as depicted in Fig. [31 

The pseudocritical points Xt° an< ^ X/T a ^ so obey scaling 
laws. Their behavior for T ~ T c and h ~ is described 
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FIG. 4: (Color online) Peak heights of the chiral susceptibili- 
ties x^r an< A X V h plotted versus the reduced mass h on log-log 
scale, for the MT1 model. The lines are linear fits through 
the calculated points which are used to obtain the critical 
exponents according to Eqs. (fT8|l - ((19}. 

by 

X p T c cc h~ x+x / 5 , (18) 
X p h c oc hM 1 -^ . (19) 

Following [iij we use these scaling relations to obtain /3, 
5, and T c for the all models and parameter sets given in 
Tab. U as illustrated in Fig. g] 

VI. RESULTS AND DISCUSSION 

As already mentioned above, all interactions yield a 
second order phase transition with mean-field critical ex- 
ponents, regardless of the strength of the interaction, its 
range, or its pointwise behavior in any particular momen- 
tum range. In particular, the presence of the 5-function 
term does not make a difference in this respect. An- 
other general observation is the fact that the function 
D(p 2 ,u)k) in Eq. @ is identically zero for all interac- 
tions in rainbow truncation. 

The values obtained for the chiral transition tempera- 
ture T c , however, are rather different among the various 
interactions. In the following we will argue how this fea- 
ture can be exploited to discern various forms within a 
given truncation. The results for T c are summarized in 
Tab. HI1 for the four interactions on a range of model pa- 
rameters well-used in meson phenomenology. 

The MN interaction has to be treated somewhat sep- 
arately, since the only free parameter in this case is D 
whose choice completely determines Q . In terms of plain 
numbers it is interesting to see that MN yields the high- 
est value for T c , followed by MR. AWW and MT give the 
smallest numbers. 

For a given interaction, T c increases with w, which is 
illustrated in panel (a) of Fig.[SJ A straight-forward inter- 
pretation of this observation goes back to uj representing 



TABLE II: Results for the transition temperature T c rounded 
to MeV for the interactions and parameter sets defined in 
Tab.Q] 



w [GeV] 


N/A 


0.3 


0.4 


0.5 


Model 


MN 








T c 


169 








Model 




AWW1 


AWW2 


AWW3 


T c 




82 


94 


101 


Model 




MT1 


MT2 


MT3 


T c 




82 


94 


96 


Model 




MR1 


MR2 


MR3 


T c 




120 


133 


144 



an inverse effective range of the interaction. In this pic- 
ture, T c is higher for an interaction, for which this range is 
smaller, i.e. the major strength in the interaction comes 
from a region defined by the larger momentum scale uj. 
However, this effect is obviously minor compared to the 
differences with regard to the form of the interaction, in 
particular the appearance of the 5-function term. 

To further illustrate the details of the differences be- 
tween the various forms for the interaction as well as to 
highlight their characteristic features, we plot T c as a 
function of three more specific properties of each model, 
namely uj D in panel (b), D in panel (c), and the chiral 
condensate in panel (d). While the attempt to quantita- 
tively correlate T c globally (i.e. across the different inter- 
actions) to any of these fails, a clear structure is visible. 
Throughout all panels of Fig. [5] the AWW and MT re- 
sults are rather close together, but clearly separate from 
MR, which again is clearly separate from MN. 

Both uj D and D investigated in the two middle panels 
have been interpreted as an "integrated strength" of the 
interaction. In fact, the latter is proportional to the inte- 
gral over d 4 q of the interaction at zero temperature for all 
interactions investigated here, if one leaves out the UV 
part Tjjy in Eqs. (IT^I) and (flUf . Fig.[S]in panel (c) clearly 
shows that there is no overall simple dependence of T c on 
D. In the MT interaction, a constant value for to D leads 
to (almost) unchanged masses and decay constants for 
ground-state pseudoscalar and vector mesons, for which 
reason uj D can be termed "integrated strength" in this 
case instead of D. Again, Fig. [3] in panel (b) indicates 
no simple dependence of T c on ujD. In both cases, the 
absence of a simple relation is exemplified already by the 
fact that three points with the same uj D, the character- 
istic feature for MT, in Fig. [5] (b) and three points with 
the same D, the characteristic feature for MR, in Fig. [5] 
(c) each correspond to three different values of T c . 

Finally, in Fig. [5] (d) we plot T c as a function of an 
order parameter of chiral symmetry breaking, the chiral 
condensate, evaluated at T = 0. In a very simple pic- 
ture, one could argue that, given a certain form of the 
dependence of the order parameter on the temperature 
(see, e.g. Fig. [3]), an enlargement of the value of the order 
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FIG. 5: (Color online) From left to right, the four panels show the results for T c as a function of: cj, leftmost panel (a); 
ujD, left-middle panel (b); D, right-middle panel (c); and the chiral condensate, rightmost panel (d). The different colors and 
symbols on the lines denote the corresponding interactions, namely: black (circles) - AWW; red (X) - MT; and blue (diamonds) 
- MR. The result for MN (single green square) is only shown in panels (c) and (d), since it does not have a value of ui associated 
with it. The insert in panel (d) enlarges the region of interest for the ^-dependent interactions. 



parameter at T = would naturally lead to an increase 
in T c . However, panel (d) clearly shows that the interac- 
tions used in our setup are not that simple. 

Overall, these observations lead to the important pos- 
sibility to establish distinct ranges of values for T c acces- 
sible for each type of interaction, which has the potential 
to rule out certain forms in a given truncation. We note 
here that various effects need to be taken into account 
beyond the simple setup presently used; we will attempt 
to quantify some of them below. At this point, how- 
ever, we can still try to attribute the clear differences 
apparent from all panels in Fig. [5] to the structure of the 
interaction. The object of interest in this respect is the 
appearance and relative strength of the <5-function term 
in g, Eqs. {JO} - (T3J): In the AWW and MT interac- 
tions, Eqs. (jTTJ) and (fT2|) . there is no such term. In the 
MR interaction, Eq. (fT3|) . it carries half the strength of 
the coupling in the sense of the integral J d 4 qQ ~ D de- 
scribed above. In the MN interaction, Eq. (fT0|) . it is the 
only term and obviously represents all of the interaction's 
strength. From this, one can say that a less pronounced 
(5-function term in the interaction will lead to lower val- 
ues of T c . 



VII. BEYOND RAINBOW TRUNCATION 

While the rainbow truncation of the gap equation pro- 
vides a simple setup for computations, one has to keep in 
mind that the necessary corrections in any given circum- 
stance may change its results considerably. Therefore 
we discuss here some of the most important effects from 
which a change, in particular to the critical temperature, 
is to be expected. 



From a phenomenological point of view one can try to 
quantify the contribution of the rainbow-truncated gap 
equation to the full (untruncated) result. Such an at- 
tempt has been made for meson properties at zero tem- 
perature in [43j and has been extended to the baryon 
sector and exemplified for several hadronic observables 
in [45] . The idea is to change the values of the available 
model parameters such that the resulting hadron prop- 
erties, e.g. the rho-meson mass, are deliberately overes- 
timated by the rainbow- ladder truncation result. In this 
way one expects corrections beyond this truncation to 
bring it to the experimental value. 

While we do not want to discuss the particular as- 
sumptions made in [43|, |4|| , the effect described there is 
relevant for our present study, since we start with model 
parameters fixed to meson phenomenology at zero tem- 
perature. To get a first estimate of the effects from a 
change of parameters at T — on T Cl we simply adapt 
the two parameter sets referred to as "A" (which corre- 
sponds to our MT2 parameter set) and "B" (which has 
increased strength in the interaction to overestimate the 
rho-meson mass as discussed above) in Ref. [45| and add 
T c to the Tab. I given there. In this table, a systematic 
ratio of ~ 0.74 is found for value in A divided by value 
in B for several hadron properties (e.g. m p , / T , mjv, ttia, 
etc.). Our values for T c in this respect are 94 MeV for 
set A and 129 MeV for set B, which produces a ratio of 
0.73 and thus fits perfectly into the picture. 

For the present discussion here the relevant point is 
that parameter changes inferred from the phenomenolog- 
ical estimation of corrections expected beyond rainbow- 
ladder truncation are in agreement with investigations at 
zero temperature. In particular, one obtains an increase 
of T c by about one third, depending on the details of the 
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corrections assumed. 

A different investigation by D. Blaschke, D. Horvatic, 
and D. Klabucar regarding the role of the Polyakov 
loop in an analogous, separable setup of the quark DSE 
at finite temperature indicates that the inclusion of 
a Polyakov-loop potential shifts T c up by ~ 70 MeV 
[46j . This is also observed in other approaches such 
as Polyakov-quark-meson studies or Polyakov-Nambu- 
Jona-Lasinio-type models, see e.g. [47} and references 
therein. 

Finally, we briefly quantify corrections from inclusion 
of more structure in the quark self-energy, in particular 
in the dressed quark-gluon vertex. Ideally, one would 
together with the gap equation also solve the Dyson- 
Schwinger equations for the gluon propagator and the 
quark-gluon vertex — the two unknown ingredients in the 
DSE for the quark propagator — as self-consistently as 
possible at any given finite temperature. Such investi- 
gations have already been started at T = (see e.g. [48[ 
and references therein) and represent a natural and very 
promising, but at the same time rather involved path for 
future extensions of the present study. 

In a more immediate way, one can identify corrections 
to rainbow truncation from the terms in the quark-gluon- 
vertex DSE and attempt to incorporate these into the gap 
equation. Pion-loop contributions have been investigated 
at T = in [49| . Ref. [5(| reports preliminary results of 
an analogous study at finite temperature, which indicate 
a reduction of T c of about 5% upon the inclusion of pion- 
loop, but not scalar-meson-loop corrections in the quark 
self energy and no change of the (mean field) critical ex- 
ponents. 

Another effect, coming from the combination of two 
other vertex corrections, can be estimated by an exten- 
sion of a simple MN-study of quark and meson properties 
presented in [5l| to finite temperature. This is a straight- 
forward procedure, which has many aspects and will be 
investigated in full detail elsewhere (52j . For the moment 
we note that a first estimate of shifts in T c are of the 
order of 15%. Another interesting observation is that in 
the course of dressing the bare quark-gluon vertex one 
can arrive at results with the possibility for solutions of 
the gap equation with D(p 2 ,ujk) ^ 0. 

To conclude this section, we remark that corrections 
to the rainbow truncation of the quark DSE should not 
be viewed as corrections to a mean-field approximation; 
corrections in the quark self energy in our present cal- 
culations happen at the level of the structure of the 
quark-gluon vertex. We would also like to emphasize 
that, although not demonstrated explicitly here, for fi- 
nite current-quark masses we obtain a smooth crossover 
instead of a second-order phase transition. 



VIII. CONCLUSIONS AND OUTLOOK 

We have investigated the impact of various forms 
of an effective quark-gluon interaction in the rainbow- 
truncated quark DSE in QCD at finite temperature on 
the critical temperature of the chiral phase transition. 
The established mean-field behavior of this transition in 
the present truncation is confirmed throughout. Regard- 
ing non-universal aspects of the transition, it is apparent 
that there is no simple overall relationship between the 
distribution of the infrared strength of the interaction 
and the critical temperature, although 1) a trend can be 
extracted identifying the part of the interaction strength 
carried by a (5-function term in some of the models as 
a qualitatively distincitve parameter, and 2) within the 
various models simple relations of T c and relevant model 
parameters do exist. 

Changes of T c due to various corrections to the present 
truncation are quantified and put into context. In addi- 
tion to effects already described in the literature we esti- 
mate particular types of corrections in the present work 
using phenomenological arguments on one hand as well 
as a simple model-extension of the quark-gluon vertex on 
the other hand. Alltogether, the various corrections put 
our results in a consistent perspective with regard to the 
results for T c obtained in other approaches. 

Furthermore, the present results clearly show that the 
chiral transition temperature is an excellent candidate for 
detailed studies of various aspects of the strong interac- 
tion. In particular it distinctively relates to both assump- 
tions about as well as particular and detailed properties 
of the effective interaction between quarks and gluons in 
the DSE setup. Further studies beyond the present inves- 
tigation will shed more light on this relationship as well 
as on how the mean-field behavior of the chiral phase 
transition may change beyond the current truncation. 
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